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ABSTRACT 


This  note  examines  the  power  required  for  heave  attenuation  of  SEV,  using 
either  cushion  air  dumping  or  flow  modulating  fans.  Shortage  of  time  has 
prevented  a  complete  analysis  of  some  aspects,  particularly  for  the  modulated 
flow  case  which  is  of  most  practical  interest.  Power  requirements  have  been 
calculated  for  100%  heave  attenuation  only  in  the  latter  case  and  the  values 
obtained  are  naturally  very  large.  More  work  needs  to  be  done  on  the  partial 
alleviation  case,  which-is  of  more  practical  interest. 

The  paper  first  establishes  the  cushion  volume  change  in  a  sinusoidal  sea 
and  from  this  calculates  the  instantaneous  pressure.  The  so-called  "compresi- 
bility  terms"  increase  in  size  by  a  factor  I (where  L  is  the  cushion  length) 
when  Froude  scaling,  but  this  is  due  to  speed,  not  size.  They  are  also 
dependent  on  the  fan  characteristics,  and  vanish  if  3p/3Q  ■  0.  They  do  not 
exist,  therefore,  in  the  case  of  100%  heave  attenuation  by  fan  modulation. 

Even  when  they  do  exist,  it  is  not  clear  that  they  necessarily  increase 
the  dynamic  cushion  pressure  excursions. 

Partly  for  this  reason,  and  partly  because  compressibility  terms  complicate  the 
analysis,  the  heave  motion  is  that  solved  for  the  incompressible  case  only. 
Simple  equations  are  obtained  for  the  power  dissipated  in  heave  attentuation . 

The  power  requirement  of  a  fan  flow  modulation  Heave  Attenuation  System  (HAS) 
appears  very  sensitive  to  the  total  head  losses  which  occur  between  the  fan 
and  the  cushion,  and  to  vehicle  speed. 

It  is  believed  that  this  note  establishes  the  feasibility  of  studying  HAS 
dynamics  analytically  and  thus  establishing  broad  trends  and  identifying 
potentially  optimum  solutions.  But  there  have  been  insufficient  time  available 
to  go  all  the  way  to  these  final  objectives  in  this  first  look  at  the  problem. 


INTRODUCTION 


The  purpose  of  this  paper  is  to  investigate  the  approximate  air  supply 
requirements  of  an  ACV/SES  in  regular  waves,  in  a  simple  yet  realistic  manner. 
One  of  the  first  such  analyses  was  by  Beardsley1,  who  considered  only  the 
"platforming"  case  over  sinusoidal  waves.  Subsequent  analyses  have  included 
many  variables  so  that  solutions  can  only  be  obtained  numerically.  It  is 
the  intent  of  this  present  paper  to  bridge  the  gap,  to  a  certain  extent, 
between  the  rather  extreme  simplifications  of  Beardsley  and  the  more  sophisti¬ 
cated  computer  studies,  without  losing  the  virtues  of  simplicity  and  closed 
form  solutions. 


Since  our  purpose  is  to  obtain  approximate  "order  of  magnitude"  results,  we 
first  note  (Figure  1)  that  the  parameter 


or 


P0Q0  - 


So 

s 


constant,  approximately 


Here  p  and  0  are  the  equilibrium  values  of  the  cushion  pressure  and  air 
volume°flow.  As  shown  in  Figure  1,  most  vehicles  fall  in  the  range 

p  0  . 

4  *  "T2,  >  T  <ft  lb/sec) 

p  0  is  the  measure  of  the  energy  lost  in  the  leaking  cushion  air.  Thus, 
awell-sealed  SES  tends  to  have  low  values;  ACV's  with  a  large  clearance  gap 
have  high  values. 


Figure  2.  Basic  Geometry. 


We  assume  a  craft,  at  forward  speed  u  ,  heading  into  a  sinusoidal  sea.  The  craft, 
planform  is  rectangular,  of  length  L  and  beam  B,  and  is  assumed  to  be  incapable 
ot  pitching.  The  "air  gap"  is  negligible  in  relation  to  the  cushion  height  H.  5 
The  local  wave  elevation  is 


h 


h  sin  2v 
w 


[f-^0 


(1) 


At  a  given  time  t,  the  cushion  volume  will  be 

L/A 


V  =  LBH 


/r  (u  +U  )  "1 

sin  2ir  [  £  +  °x  —  tj  d(x/A) 
o 

h  <vu»’  .1  .  <vu«>  A 

Lx  *  — x—  *J  - cos  2*  — r—  ‘I 


(2) 


h  BA 
w 

2ir 

h.  BA 


i  cos  2ir 


=  LBH  — /  sill2  2ir(L/A)  +  [cos  2ir(L/A)  -  l]2  sin  [ -2-^Q-  t  ♦ 


(where  is  a  phaso  angle,  and  is  not  important,  and  U  =  u  ♦  u  .  the  sneed 
relative  to  the  waves.)  Continuing  the  redution:  o  w  * 


V 


(3) 


h  BA  _ 

=  LBH  -  -~-y/2[\  -  cos  2tt(L/A)]  sin[2ir  (U/A)  t  +  ♦j] 

We  pause  here  to  note  that,  for  comparison  with  simple  "piston  theory,"  the 
displacement  (5)  of  the  "wave  piston"  is  given  by 


£  '  sl"[2*W»>  *  *  «i] 

-  sin^-WX)  t  .  ♦,]  (4) 

2 

This  is  the  result  given  (for  example)  by  Mantle. 

Returning  to  equation  (3) ,  the  rate  of  change  of  cushion  volume  is 

3F  •  LB  j§f  '  hw  (u/i)  “S[2'(UA)  1  *  *l]  j 

-  Lb[^  -  nhwF(L/A)  COS  (fit  +  fj)]  (5) 

when  fi  =  2ir(U/A),  the  frequency  of  wave  encounter  in  rads/sec. 

“d  pw‘)  ■ 

We  shall  show  later  that  the  heave  term  dH/dt  is  generally  negligible  by 
comparison  with  the  cos  fit  term.  The  phase  angle  may  usually  be  ommitted, 
since  it  can  be  made  zero  by  a  suitable  selection  of  zero  time. 


111)  IJ»«j 


PC  +  P- 
P0  +  P- 


QF 

Qo 


THE  CUSHION  PRESSURE 
absolute  cushion  pressure 

absolute  cushion  pressure  under  equilibrium  conditions 
fan  volume  flow  rate 

fan  volume  flow  rate  under  equilibrium  conditions 

mass  of  air  in  the  cushion  under  equilibrium  conditions 

air  mass  flow  rate  into  the  cushion  (sum  of  flows  in  and  out) 

cushion  volume  under  equilibrium  conditions 

equilibrium  cushion  density  =  ®0/vQ 

o  at  an  equlibrium  conditions 


l 


At  any  time  t  the  density  of  the  cushion  air  is 


For  adiabatic  conditions 


V  ^ 

■  T  <  1  *  dt> 
v  I  Pc  *  p.\Vt  , 

Vo'Po  +  P  J 


(6) 


(7) 


Differentiating  with  respect  to  time 

_m_  _  J_  dV  /Pc  +  P»V/Y.  V  /Pc  +  P-^-1  /  1  \dPc 

mo  "  Vo  dt  VP0  +  P J  +  YV0'P0  +  P •'  \P0  +  pJ  <*t  (8} 

Since  m  *  p  V 


1/y 

m  =  po  /  Pc*P»  \  fdV  (  V  dp  ~1 

p»  p»  '  P0+P»  '  Ldt  Y(PC+PJ  dt  J 
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=  (Pg  ^  l\1/Yfdv  -  v 

VP»  1)  Ldt  Y(P  +pj  dtj 


(9) 


dV 


for  incompressible  flow  (y  -*■  ®) 

If  it  is  assumed  that  p =  p  +  Ap  and  that  Ap/p  «  1,  then  equation  (9) 
becomes  c  o 


i  -  £  ♦!.) 

P„  \Pco  / 


dV  +  V  d£ 

dt  *(P0  +  PJ  dt 


(9a) 


The  flow  out  of  the  cushion  through  the  leakage  area  (a)  is 


Qj  ■  a/(2/p)  pc  *  a/(2/p  J  (po+Ap) 


The  flow  into  the  cushion  from  the  fan,  (neglecting  inertia  terms)  is 

Of  ■ 


Making  these  substitutions  in  equation  (9) 

%  +  Ip  AP  -  a/(2/Pj(P0+AP)  =  (?T  +1) 


1/Y1 


dv  .  _ 

eft  +  y 


V  (Po+P»)1/Y_1  d 

"  i7y 


& 


(10) 


here  V  and  dV/dt  are  given  by  equations  (3)  and  (5) .  This  equation  cannot  be 
solved  explicitly  for  Ap,  even  though  it  is  first  order,  because  of  the  square 
root.  An  approximate  solution  is  discussed  in  Appendix  II. 

Some  convenient  non-dimensionalizations  are 


%  3£  _1_  dV 

P0  3Q  *  ^  =  Qodt 


(10a) 


*  3Q/3p  is  not  a  true  partial  derivative  (as  witness  its  subsequent  inversion) 
but  the  notation  used  is  the  most  familiar  in  the  literature.  In  retrospect 
thd  more  gengral.  expression' 

AP  =  cQ  +  CjQp  +  ;c2QfZ 

-  might  lead  to  a  more  general,  and  thus  more  useful  result. 
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Equation  (10)  then  has  the  alternative  form 


♦  E 

?Po 


/I  +  Ap/p. 


=  C(l+PD/pJ 


1/Y 


U+P0/P«,) 


l/Y-1 


V  dp 

YQo  dt 


(11) 


The  terms  in  the  brackets  are  the  so-called  ''non-scaling*'  and  "compressibility" 
terms.  (For  example,  see  Lavis,  et  al3).  The  curly  bracket  is  equal  to  zero 
for  incompressible  flow  (y  -*■  <*»),  the  "non-scaling**  term  equals  unity. 

If  fan  modulation  is  employed  to  eliminate  heave,  dp/dt  =  Q,  then  an  incompres¬ 
sible  flow  analysis  is  in  error  by  the  factor  (Pq/P,,,  +  l)1***  !•©• 


if  p  10  SO  100  200  lb/ft2 

o 

(p0/pw+D1/Y  1.003  1.017  1.034  1.067 

This  is  small  enough  to  be  neglected  in  an  analysis  of  the  present  type.  For¬ 
mally,  one  can  account  for  it  by  assuming  a  wave  height  which  is  lower  by  this 
factor. 


On  the  other  hand,  without  attenuation, 
for  incompressible  flow,  approximately 


so  that 


we 


shall  find  (equations  18  and  18a)  that 


Equation  (10)  then  becomes,  approximately 


!  +  i^E  _  /f*Ap/p0 
s  *o 


<>*p0/pj1/’r 


jl  ds  ! 

Vi+p/pJ\2-«/  yqd  at  J 


Substituting  equations  (3)  and  (5)  for  V  and  dV/dt,  and  assuming  that  craft 
heave  motion  is  small  compared  to  wave  motion,  so  that  LB(dH/dt) 
is  negligible  compared  with  dV/dt 


1  •  if  -  ■  '''"v*.]  “i- *  •  “■ 01 

41 


-  F (L/X)  sin2ftt 


(12) 


where 


fi  y  h 

6  =  -q^FCL/A)  «  2,^-LB-f 


=  2irU 


-f— ^ 

X  'Vo/ 


Several  points  are  apparent  from  this  result: 

1.  The  "compressibility  term"  gives  a  phase  shift  and,  even 
on  these  simple  assumptions,  introduces  both  a  constant 
term  and  a  second  harmonic. 

2.  Since  G  is  independent  of  scale,  there  is  no  significant 
scale  effect  for  constant  V0/Vm> 


but 


3.  Since  G  varies  with  absolute  speed,  it  does  become 
important  when  Froude  scaling;  it  increases  as  /l,  so 
comparing  a  30th  scale  model  to  full  size,  G  is  five  times 
as  much  in  the  full  scale  case.  And  since,  in  Froude 
Scaling,  cushion-pressure  must  vary  with  length,  the  total 
variation  is  as  l$/2,  for  a  total  factor  of  164  in  the 
example  of  a  l/30th  scale  model. 

4.  The  "compressibility  term"  will  not  be  important  if  the  fan 
characteristic  parameter  (  is  small.  It  will  be  zero  if 
3p/3Q  =  0.  In  practical  applications,  the  parameter  C  varies 
between  0  and  -2. 

5.  Because  of  complex  phase  shift  effects,  the  compressibility 
terms  may  either  attenuate  or  increase  the  heave  acceleration, 
depending  on  the  precise  values  of  the  parameters.  We  show 
in  Appendix  II  that  if  the  heave  velocity  dH/dt  can  be  con¬ 
sidered  negligible,  compressibility  will  always  attenuate  the 
ride.  It's  therefore  important  to  understand  this  effect,  and 
design  (so  far  as  other  factors  allow)  for  a  favorable  effect. 

6.  The  phase  shift  in  peak  pressure  will  generally  result  in 

an  increased  heave  attenuation  system  (HAS)  power  requirement. 

In  view  of  these  considerations,  the  work  which  follows  will  ignore  the  "com¬ 
pressibility"  term,  and  consider  only  incompressible  flow. 
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HEAVE  MOTION  FOR  INCOMPRESSIBLE  FLOW 


For  incompressible  flow  equation  (12)  becomes 


V  p  (Po  +  Ap)  =  %  *  3?  Ap 


3p  r  dt 


Again,  let 


5  po  3Q 


J_  dV 
Qo  dt 


(i.it)  . 
p  M  Po; 


And  since 


/I  ♦  4e  =  i  ♦  I^E- 


po  «  po 


Squaring  both  sides  we  obtain  the  quadratic 


(£)*  * [2tl  - »  -  «  (&)• C1  -  »2  - 1  •  • 


The  plus  root  is  taken  in  this  case. 


{('  -  f)  -  *] 


1  -  f  ♦ « 


i 
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PEAK  ACCELERATION  IN  Qt  OR 


Figure  3.  Exact  and  Approximate  Solutions  to  the  Equation  for  Cushion 

Pressure.  A p/pQ  =  f(5)  so  Ap/pQ  =  f(5Bax)  (Equations  18  and  18a) 
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If  we  formally  restrict  the  analysis  to  Ap  «  p 


2cKh 

[2  -  C] 


2^p  nK 
o  w 

[2  -  d 


?a  Js-  =  g  _ 
Qo  2  -  «  g  p, 


w_  /  2C  \ 

A.  I2  -  V 


Zlt  [- Zxt 

Ap  «  e  I  /e  sin  Ot  dt  +  C 

z  zt 

=  -  —  cos  at  +  Ap  e 
a  o 


Thus  for  steady-state  conditions  (no  transients) 


^  '  •(x')  p(v)M%) cos  “ 


So  finally,  noting  (19)  and  comparing  equation  (25)  with  the  equivalent 
"piston  motion"  of  equation  (4)  we  get  the  rather  simple  result 


—  s  craft  acceleration  __  / gS  \  2g 
6  B  piston  acceleration  “  yQQ^  [2  -  t] 


Comparing  (25)  with  (18a)  we  see  that  they  are  identical,  if  the  dH/dt 
term  in  equation  (5)  is  dropped.  So,  to  a  first  order,  we  may  neglect  the 
heave  motion  of  the  vehicle  in  computing  the  air  flow  into  and  out  of  the 
cushion. 


*  From  (18a)  and  (5) 


<*)  ■  i  s-  ¥  •MW— 
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The  transient  term  in  equation  (24)  is  independent  of  the  forcing  term  Z2» 

and  if  h  =  0,  Z  =  0  and 
w  ’  2 


Ap  =  ApQe  1  (27) 


where 


•(I) 


2C 


[2  -  C] 


Thus  the  motion  will  be  unstable  if  0  <  c  <  +2;  a  known  result.*  For  all  other 
values  Z^  <  0  and  the  transient  decays,  so  that  the  motion  is  stable. 


Heave  Attenuation  by  Dumping 

Let  6  »  IK  and  the  suffixes  o  and  D  refer  to  dump  valve  shut  and  open. 
Ap,  is  the  pressure  differential  when  the  dump  valve  opens  (at  fit  «=  9j) 
and  Ap2  the  valve  when  it  closes,  as  indicated  in  Figure  4. 

By  definition 

•>0  -  “TF^T  [Equation  (16)] 

Also 

%  "  &D>/p  po 


where  aD 


is  the  normal  leakage  area  plus  the  dump  valve  area. 


[corresponding  to  Z.  in 
equation  (23)] 


*  A  result  first  noticed  by  Walker 
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CUSHION  PRE8SURE 


Sh  (2 

Z  *  * 

0  % 


Sh  0 


ZD  - 


€  ^ 


[corresponding  to 
Z2/P0n  in  (28) 

equation  (23)] 


when  the  dump  valve  is  closed 


AP  *  -P0Z0  cos  6  ♦  Ap2e 


zo^*2 


and  open 


AP  =  -Pqzd  cos  0  +  Ap^e 


zDAtl 


e  -  e. 


At 


2  a 
6-6, 


At 


1 


So  at  the  point  of  opening  the  dump  valve,  at  0j 

AP^  a  -P0Z0  cos  6 j  +  Ap2  Jzo^n(0i"®2)] 
at  the  point  of  closing  the  dump  valve,  at  6^ 

Ap2  “  “P0Zd  cos  e2  +  Apje^D^  ^  -0l)J 

By  substituting  in  (30a)  for  Ap2  from  (30b) 


[corresponding  to 
equation  (24)] 


(29) 


(30a) 


APi 


-p,Z  cos  0 
*0  o 


[  b 

1  +  -P0ZD  cos  02  +  Aple 

m 


TT<VV 


irW 


(30b) 


/ 


r  z  /z„ 

\ 

Z 

APi 

S.  L2.  .  3 

1  -  e°  'Z0 

*  -Po 

IT  <Ve2> 

ZQ  cos  0j  +  Zp  cos 

(31) 


Similarly,  substituting  for  Apj  in  (30b) 


Ap2 

.  -  _  -  7 

.  „  i  (W1 

cos  0j  +  cos  @2  e  1 

Po  "  ° 

.  t  (*-£)  <VV  . , 

(32) 


So  by  selecting  values  for  0.. ,  6-  and  Zn/Z  we  can  compute  Ap.  and  Ap« , 
and  hence  evaluate  equations1 (29J.  BecausS  of  the  exponential  terms 
we  cannot  easily  obtain  maxima  analytically.  The  portion  of  the  solution 
in  which  Ap  will  occur  is,  from  (29)  and  (31) 


A£ 

Po 


r 

ZD  ■£TCei“02) 

V  ) 

=^-  COS  0  ♦ 
Lo 

< 

cos  +  cos  ®2  6 

o 

eir(6-vj 

( 

(33) 


It's  clear  that  the  maxima  must  be  sought  graphically  or  iteratively;  which  ii 
easy  to  do  by  computer,  of  course.  Note  that  the  variables  are  four  in  number 


/ 


i.e.  &.  =  f  (0J  82  zo  zD) 

*o 

A  Solution  in  Which  the  Transient  Terms  are  Ignored 


In  this  case,  when  the  dump  valve  is  open  [from  equation  (29)] 


z^  cos  e 


(®j  <  0  <  e2) 


(34) 


/*E\ 

\Po/ 


-  -2, 


max 


And  the  attenuation  ratio  is 


(aPmax)D 

fPmax)T 


We  open  the  dump  valve  when  Ap  is  equal  to  the  maximum  value  of  Ap  which 
will  occur  with  the  dump  valve  open. 


i.e.  /AR\ 

'  °'max 


APi 


ZQ  cos  02  =  -ZD 


cos  e. 


sin  8, 


-Z, 


Also,  from  Figure  4 


e2  s  Zir  -  0! 


so  that 


sin  02  ■  -sin  0j 


The  instantaneous  power  lost  through  the  dump  valve  is 


PD  ■  ^  “D  f (Po  *  ap)  ‘  15  0(aD  -  *>(f  ^  (‘  *  1^) 

■  *sp  (ap  -  a)||-  P0j  |l  ♦  J  Zp  cos  0  j  [from  equation  34)  (35) 
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Averaged  over  one  cycle  this  is,  as  a  ratio  of  equilibrium  power 


DAV 

P 

o 


'  (t  -  *)  5T  l'e2  -  V  *  I  Vs1"  °2  -  sin  61 


)] 


-  (t  -  *)  h  V1  -(rf 

I-  o  J 


(36) 


Thus  the  power  ratio  is  not  a  simple  function  of  area  ratio,  ap/aQ»  but  depends 
upon  the  attenuation  ratio: 


and 


ZD  /  \  * 

T~  [-  peak  acceleration  with  HAS  j 

o  \  peak  acceleration  without  HAS/ 


/_  leakage  and  dump  area} 
a  \  leakage  area  / 


-  ■* '  4  *  ■«**  ■ 


There  is  therefore  no  simple  way  to  present  this  result  unless  we  assume 
CD  =  V  S0  that 


% 

a 


(37) 


Then 


PDAV  .  V  +  V 
V  To’  To’ 


(38) 


•HAS  n  Heave  Attenuation  System. 


/ 


-  and  only  the  second  term  depends  upon  the  details  of  the  ship  and  waves. 
These  two  terms  are  plotted  in  Figure  5.  The  second  term  may  generally  be 
neglected,  except  under  very  rough  conditions,  when  large  accelerations 
are  experienced  with  the  HAS  activated. 


Attenuation  with  a  Variable  Supply  Fan 


Heave  attenuation  can  also  be  provided  by  a  variable  pitch  fan,  or  any  one  of 
several  other  ways  of  modulating  the  fan's  flow.  In  this  section,  we  assume 
that  the  modulation  is  loss-free. 


The  flow  required  from  the  fan  is 


Q  =  ^  +  a^"( Pd+ap)  (a  vers*on  equation  13) 


JF-'W1 


4 


*  *r 


The  total  head  rise  through  the  fan  itself  is 


AH  *  v 


(l  ♦  5E.)*  (1  -  nj.5  %  p(£)2 


(39) 


(40) 


where  (1  -  nj)  is  the  total  head  loss  divided  by  the  dynamic  head. 
Let  np  =  the  fan  efficiency  (=0.8) 


nL  »  the  total  lift  system  efficiency  (=  P0Q0/P  “  0.4) 


[p„  *  (»  -  "t)*  "(t)  ]('  - 


(from  equation  40) 


0 


02 


0.4 


1.0 


Figure  5 


i 


HEAVE 


ATTENUATION  ■ 


0.6  0i8 

PEAK  ACCEL.  WITH  H.A.S. 

PEAK  ACCEL.  WITHOUT  H.A.S. 


Power  Required  for  Heave  Attenuation  by  Dumping. 

(The  tj  parameter  must  be  multiplied  by  G  before  substitution 
in  the  equation  for  AP/P  .)  " 


(41) 


/ 


«P0  *ay 


Reverting  to  equation  (40)  the  instantaneous  power  required  is  given  by 


HpP  ■ 


“» ■  [*(-£)*  *'•(£)*] 


Let  5  =  ^ 


1_  dV 

%  dt 


as  before 


(42) 


*  1+^2. 


(42a) 


Then 


JSL 

% 

_npP 

^0% 


(p2") 

'  o  /  max 


«  (€  +  q) 


(from  equation  39) 


=  [q  ♦  §(£+  qHU  +  q) 
=  q(C  +  q)  +  §U  +  q)3 


(43) 


This  implies  that  when  AH  >  0  and  Q  <  0  (or  visa-versa)  the  fan  abstracts 
power  from  the  air  like  a  windmill,  which  may  not  be  realistic  unless  a  flywheel 
is  used.  Note  from  (39)  and  (40)  that  although  Q  can  be  negative,  AH 
cannot.  Thus  the  equation  (40)  relationship  for  AH  requires  modification 
when  Q  0. 

Now  for  sinusoidal  waves  we  have,  from  (20) 


C 


J_  dV 
% 


K  cos  ftt 
w 


Kw  cos  6 


(44) 


and 


(48) 


1 


and 


ip 

HAS  _  Power  to  Eliminate  Heave 
PQ  -  Equilibrium  Power 


(49) 


Note  that  if  n,  =  rip.  no  power  would  be  required.*  The  more  efficient  the 
lift  system  duct  ana  diffusion  system,  the  less  the  penalty. 

Without  heave  attenuation  the  maximum  acceleration  would  be  (from  25) 


A*max 

Po 


A 

P  =  Zo 


b  K*f(0  say 


(50) 


So  if  z  is  the  heave  attenuation  ratio 
P  •  *Kf(C) 

making  this  substitution  in  (47)  putting  §  ■  (hp/n,  -  1)  (equation  40)  and 
rearranging  gives 

OpPaV  Op  3  /Op  \  2  2  2 

p£qT  -  ^ +  2  (^  -  V  V  -  Ifwr  -  %  Kw2zf(C)[l  -  *f(0]  (51) 

*  Strictly  true  if  rip  ■  1.0  or  if  there  is  no  reverse  flow  through  the  fan.  When 
rip  <  l.£  and  reverse  flow  occurs,  the  fan  absorbs  too  much  energy,  by  a  factor 

of  1/rip  in  this  analysis. 


/ 


As  nFPav^po^o  =  e<luilibrium  value  (=  np/nL)  plus  the  increment  ^h/VS^o^o  »  we 
see  that 

^  "k2  |3(‘  -  ^)  1'  -  2£C‘»2  -  ^  0  -  l«0)  j  tS2) 

This  is  finite  even  for  no  attenuation  (z  =  1) . 

If  z  =  1  and  nL  =  np 


-  h  Kw2  [1  -  f(0]  (52a) 

-  implying  a  reduction  of  poyer  in  waves  if  f(c)  <  1,  which  it  usually 
is.  But  for  practical  cases,  there  is  an  increase  in  fan  power  due  to  waves 
if  (from  equation  52) 


3  (l 


AP 


waves 


1  -  zf(0 


> 


(“ 


(53) 


-  which  in  general  will  always  be  true. 

For  the  purpose  of  obtaining  a  rough  estimate  we  consider  only  the  case  of 
complete  alleviation,  so  that,  from  equations  (44a)  and  (48),  the  additional 
power  required  is  '  ~ 


AP, 


HAS 


AP, 


HAS 


(lb  ft/sec) 


(ft/sec) 


(54) 


(55) 


since  ft  ■  2w  U/A 


This  is  plotted  in  Figure  6  for  typical  values.  It's  clear  that  complete 
suppression  of  heave  is  much  too  expensive  of  power  at  high  speed  on  the 
assumptions  used  in  this  analysis.  We  should  examine  the  implications  of 
equation  (52)  more  closely  to  see  what  the  trade-offs  are  for  partial 
alleviation, 
is  very  severe. 

For  comparison,  simple  (zero  leakage,  loss-free  dumping)  piston  theory  gives, 
for  100%  alleviation 


(ft/sec) 


(56) 


Beardsley's 


1 


wave 


pumping 


theory  gives 


2ir  hw 
\  x 


(ft/sec) 


(57) 


Both  of  these  give  much  lower  power  estimates  than  equation  (55)  for  typical 
values.  Yet  on  the  other  hand,  equation  (55)  predicts  no  power  penalty  at 
all  if  h .  ■  n  .  That  is,  for  zero  duct  and  diffusion  loss;  as  noted  earlier 
This  is  ^because  our  equations  simulate  "energy  recovery"  by  the  fan  when 
the  flow  is  negative;  also,  equations  (56)  and  (57)  do  not  allow  for  duct  losses. 


It  would  seem  that  minimizing  duct  and  diffusion  losses  are  of  paramount  impor¬ 
tance  to  maximizing  the  efficiency  of  any  fan  modulation,  heave  attenuation 
scheme. 


EXTRA  FAN  POWER  REQUIRED,  BHP/LB  OF 

GROSS  WEIGHT 


Figure  6.  Extra  Fan  Power  Required  for  Complete  Suppression  of  Heave  in 
Long  Waves,  in  Short  Waves,  Multiply  Power  by  fF (L/X) 3  . 
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APPENDIX  I 

FAN  CHARACTERISTICS 


It's  usual  to  express  the  total  pressure  rise  (AH)  across  a  fan  as 


3AH  3  AH  3<fr  pto  3<> 

3Q  *  9*  JQ  3*  *  D  3* 


and 


%  3  AH  *o  3j> 

AHo  <P0  #♦ 


We  are  interested  in  cushion  pressure  p,  which  is  related  to  AH  by 
p  *  AH  -  (1  -  Hj.)  j  p(j) 

where  n_  is  the  total  head  efficiency.*  This  equation  assumes  that  the 
plenum  (cushion)  velocity  is  negligible  compared  with  (Q/A). 

*  Head  loss  after  the  fan/Hp  (Q/A)^  where  A  is  the  fan  duct  area 


A  common  case  (and  also  the  lower  limit  on  nT)  is  for  the  diffusion  process 
to  approximate  the  Borda-Camot  ’’rapid  diffusion”  case  for  which 


a  -  v  * 


-  2  |  (1  -  A/S) 


-  where  A  is  the  total  fan  area,  and  S  the  cushion  area. 


APPENDIX  II 


THE  COMPRESSIBLE  EQUATION  FOR  Ap 


33 


If 


Jl  -  Cl  *  P0/Pj1/Y  J2  -  (P0/P.)(l  *  P0/pJ1/T'1 

P  =  Ap/po  V  *  V/yQq  C  =  (1/Q0) (dV/dt) 

Then  equation  (11)  becomes 

1  +  i  P  -  &V  =  JjC  +  VJ2  (II. 1) 

Here  £  and  V  are  the  driving  terms,  and  may  be  quite  general.  In  our  present 
case 

V  =  Vq  -Vj  sin  fit 
* 

*  dV  ^  dV 

“  dt  "  yQo  dt  =  “  YVV1  cos  01 

£  =  -  YfiVj  COS  fit 


Making  these  subsitutions  in  equation  (II. 1) 

A  • 

J2(VC  -Vj  sin  fit)  ^2.  +  .  I  p  =  I  +  JjYfiVj  cos 


fit 


(II. 2) 


Several  simplifications  immediately  come  to  mind.  For  low  accelerations  p  «  1, 
so 

/l  +  P  *  1  +  i  P 


The  same  assumption  implies  Vj  «  VQ.  Thus 

’  T^r1'05  01 


(II. 3) 


or 


3?**  =  B  cos  fit 

.*.  p  *  e“At  B  f  eAt  cos  fit  dt  +  transient  terms 
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A 
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M 

i 


i 


■4 


| 


if  e  -  ot,  t  -  e/ft 


Substituting  for  A  and  B 


The  incompressible  flow  result,  equation  (25),  may  be  written  as 

»  /  2c  \ 

Po  =  C0S  nt 

=  (ftrvp/  (I .  I)  cos  nt 


(Also  the  result  of  putting  dp/dt  =  0  and  Jj  =  1  in  equation  I 1.3). 
Thus  the  relative  amplitude 

(Ap/P0)comp 

r  -  (Ap/p  ) . 

mcomp 

r  ■*  Jj  as  c  +  0,  a  result  we  have  seen  already. 


I  ' 

l 

Li 


(II. 4) 


(II. 5) 


(II. 6) 


(II. 7) 


The  phase  angle  can  be  expressed  as 


rV  n 

cos  +2  =  ryVjfl  = 


Generally  speaking  Jj  »  1,  J2  e  p^p^.  Also 


“Vo 


£o  V0  2irU 


(II. 8) 


(II- 9) 


(II. 10) 


which  is  similar  in  form  to  G/y  in  equation  (13). 

Note  that  r  is  always  less  than  unity,  according  to  this  theory,  in  direct 
contradiction  of  earlier  results.  But  in  fact,  this  follows  from  Equation  II. 3. 
Whenever  a  first  order  linear  equation  is  driven  by  a  sinusoid,  its  amplitude 
will  always  be  less  than  the  value  obtained  when  the  differential  is  neglected. 


APPENDIX  III 


A  SIMPLE  CHECK  ON  FAN  MODULATION  POWER 
FOR  100%  ATTENUATION 


Suppose  the  cushion  pressure  is  held  constant  at  the  trim  value  po,  so  that 
the  total  head  rise  through  the  fan  must  be 


AH  =  PD  +  Cl  -  n^)  j  P(Q/A)2 

The  instantaneous  fan  power  requirement  will  therefore  be 
nFP  =  AHQ  *  p0Q  +  (1  -  nT)  \  P(Q/a)2 
Let 

Q  =  (^(1  +  Kw  cos  fit)  =  (^(1  +  Kw  cos  0) 

Now  in  equilibrium, 

"fpo  -  Pofyotv^2  v1  -y 

7»CV«2  V1  -  V  -  "FPo  -  ro% 

°  PoCnF  ■  V 

Thus  the  instantaneous  power  can  be  expressed  as 

npp  =  P0Q0(1  ♦  *w  cos  0)  ♦  P0(np  -  nL)(l  +  Kw  cos  0)3 


or 


^  (1  +  cos  0)  ♦  (!  -  nL/nF)(l  +  Kw  cos  0)3 


-  ^  (1  +  K*  cos  6)  ♦  (1  -  nL/nF) (1  ♦  3Kw  cos  0 
F 


♦  3Kw2cos20  +  Kw3  cos3©) 


Thus  the  average  over  one  cycle  is 

n, 

n 
O 

which  is  the  same  as  equation  48  in  the  main  body  of  the  report. 


Pav 

P 


^  +  I  Kw2  (1  "  W 


(111.1)  | 

(111. 2) 

(111. 3) 


1 

i 

j 

S 

i 

i 

(III.4) 


i 


(III. 5) 


4 

i 


(HI. 6)  j 


.i 

* 


i 


i 


